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Given a fixed initial state, a desired Hamiltonian, and an amount of time, we provide a complete 
characterization of the set of Hamiltonians which perform the same action as the desired Hamiltonian 
on the state of interest. An example is provided where the desired Hamiltonian lies outside of the 
control Lie algebra, but implementable Hamiltonians exist within the characterization. 



In addition to not always being feasible, having com- 
plete control over a quantum system is not always nec- 
essary. Indeed, examples of this situation are numer- 
ous, and include quantum simulators designed to simu- 
late only certain tasks [U [2] , sensors which exploit coher- 
ent dynamics to measure in the nanoscale regime O |4] , 
and most protocols in quantum communication [5H7]. If 
one has a quantum circuit and wishes to implement it, 
a quantum system over which one has complete control 
would of course work, but a quantum device capable of 
implementing only one's circuit would nonetheless suffice, 
and perhaps in a more cost effective manor. However, 
many of the tools of quantum control theory are tar- 
geted towards the cases of complete controllability, pure 
state controllability, and eigenstate controllability, which 
all emphasize a need to be able to control a given quan- 
tum system arbitrarily [5HTT] . This paper develops a new 
tool to add to the toolbox of state to state map finding 
in the case of only partial quantum controllability. 

The following is our predicament: we wish to evolve 
some state under a desired Hamiltonian which lies outside 
of our control algebra of implementable Hamiltonians for 
a specified time. The goal is to give a full characterization 
of the Hamiltonians which perform the same action as our 
desired Hamiltonian on our state of interest. Moreover, 
we will provide a physically motivated example where 
such an alternate Hamiltonian exists, and lies within our 
system's control. 

Let H = Cd denote our system's Hilbert space of di- 
mension d. Then the space of linear operators acting 
on Ji can be taken to be the d x d complex matrices 
B{H) = fll(rf) [53] • Let the matrix pi denote our fixed 
initial state of interest, the matrix Hd our desired Hamil- 
tonian, the matrix Hint the internal Hamiltonian of our 
system, and the amount of time we are interested in 
evolving for. Thus, our desired final state is given by 
Pf :— UdPiU\ where Ud — e"**"^"*. Our control algebra, 
which will be denoted by £c, is a real Lie subalgebra of 
BCH), and moreover, a Lie subalgebra of the dx d skew- 
Hermitian matrices u(d) . Since we are assuming that Hd 
lies beyond our control, we know that Cq must in fact 
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FIG. 1: The initial state evolves under Hint for a total time 
rjto, with an evolution under some member of Qn for a time 
(1 — rj)to in between. This is arranged such that the state has 
effectively evolved under Hd for a time to . 

be a proper subalgebra of u{d) . 

Fixing a fraction < rj < 1, our scheme (illustrated in 
Figure [T]) will be to evolve under the internal Hamilto- 
nian for a time ^Iq, perform control operations for a time 
(1 — 77)^0 J and then resume evolution under the internal 
Hamiltonian for the remaining amount of time ^tg. The 
purpose the the parameter rj is to facilitate the case where 
the desired Hamiltonian is similar to the internal Hamil- 
tonian, so that only a brief control sequence surrounded 
by natural evolution is necessary. An example of this 
might be removing next-nearest neighbour couplings in 
a chain of spins, desiring only the effects of nearest neigh- 
bour interaction. It can, of course, be set to zero if it is 
not needed. 

We will denote the state of our system prior to the con- 
trol sequence by p_ , and the state immediately following 
the control sequence by p^, so that 

p_ - U.ntP^Ul, (1) 
P+ - UltPfU^nt (2) 

where Uint — e^*3*°^'"'. Those Hamiltonians which 
evolve p_ into pj^ in time (1 — 77)to constitute a set which 
we will call Q-u- Our method will be to first character- 
ize the set Qjj of unitary matrices which take p_ to p+. 
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When this is done, Gn will be found by taking all loga- 
rithms of all elements of G-h . We may write 



gu = {ueu{d) I p+ = up-u^} 



gn={H e iu{d) I e 



-it„{l-r,)H 



eGu} 



(3) 

(4) 



where hl{d) and iu{d) are the d x d unitary matrices and 
Hermitian matrices respectively. 

The states p_ and p+ being unitarily equivalent, we 
may diagonalize them to the same matrix as follows. Let 
U- be a unitary matrix which diagonalizes p+ to the 
diagonal matrix D, so that p_ — U-DU^_. Then U+ 

^int^d^int^- ^ Unitary matrix which diagonalizes p+ 
to the same diagonal matrix D. 

In general, we can think of a diagonal matrix as a block 
diagonal matrix where each block is a scalar multiple of 
an identity matrix and no two adjacent blocks share their 
diagonal value in common. If T is a diagonal matrix, 
let lower-case tj be the scalar value corresponding the 
j^^ block, pj be the size of the j*"^ block, and itit be 
the number of blocks. For example, this gives us D — 
diag(d,/po)™fi. 

Without loss of generality, by swapping the columns of 
U- and [/+, we may suppose that D has been chosen so 
that di 7^ dj whenever i ^ j. Now p_ and share the 
same set of eigenvalues, namely {c?j}jl°, and with the 
same degeneracies 

We begin our characterization of Qu and Q-n with an 
easily verified matrix property. Suppose X G GIj{n) 
and that T is a diagonal matrix such that ti ^ tj 
whenever i j. Then T = XTX^^ if and only if 

X eGLipJ,...,pl^) m- 

As a result, we conclude that HUE U{d) is a uni- 
tary matrix, then D = UDW if and only if [/ G 
lA{pi , . . . ,p^^) C U{d), from which (*) in the follow- 
ing string of equivalences follow: 

W eGu ^ p+^Wp^W^ 
^ U+DUl = W{U^DUl)W'^ 
<^ D 



UIWU-DUIW^U^ 



u -.^ ulwu- 



- eU{p^ , . . . , j.„j^ 
W = U+UUl, with U e U{pf, . 



And so our characterization of Gu is just a left and right 
translation of the block diagonal group U{p^ , . . . 



by the unitaries U+ and C/1 respectively: 



Gu = U+Uipf,...,pg^)U^ 



(5) 



The next step is to characterize the set Gn of Hamil- 
tonians which carry us from p_ to p+. The difficulty in 
this step arises from the non-injectivity of the matrix ex- 
ponential function. For each U E Gu there are at least 
countably many valid matrix logarithms, and uncount- 
ably many if U has any degenerate eigenvalues. 

For any U E Gu let Vu be the set of ordered pairs 
describing the diagonalizations of J7, i.e., define 



U = VTV\ where 
V eUid-u) and 
T is diagonal 



(6) 



Since all of the T's in Vu are equal up to permutation, 
we can use the matrix property stated above in a similar 
way to derive the following more useful description of Vu 



Vt, = <^ [V'WP, P^T'P) 



P a permutation matrix 



(7) 

where {V\T') £ Vu was chosen so that t'i ^ t' j when- 
ever i ^ j. 

With all of this notation aside, we can start to talk of 
logarithms. Given some U £ Gu, a-U of the logarithms of 
U , that is, all of the matrices (skew-Hermitian or other- 
wise) which exponentiate to [/, are of the form 



A ^ VXTdieig{i(t>j + 2mkjY.^^X:^W^ 



(8) 



where {V,T) e Vu, T = diag(e*'^Oi=i: and Xt is an 
invertible matrix which commutes with T, i.e. Xt lies in 
the general linear centralizer of T, Cgl{T). It is easy to 
check that a matrix of this form exponentiates to U : 



= V XTdiag{e"*'^)X:j,^V^ = U 



(9) 



For a proof that the form in Equation [8] constitutes 
all such logarithms, the reader is redirected to Refer- 
ence [12] .We finally arrive at the following characteriza- 
tion of Gu'- 



Gn 



(1-^)^0 



V XT<i-ia,g{- 



2nkj)UXT^V^ 



UeGu, iV,T)eVu with 
T = diag(e*'^Oj=i> and 
XreCGLiT), k^eZ 



n iu{d) 



(10) 



Note that we have multiplied the form in Equation [8] by (— i(l — ri)to) ^. This is simply so that ii H E Ghi 
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then e-<^-''^*°" = U rather than = U. 

To see that the intersection with ivi{d) is necessary, 
in other words, to see that a logarithm of a uni- 
tary is not necessarily skew-Hermitian, observe that 
our characterization of logarithms tells us that A = 
(oi)(o2°i)(oi) exponentiates to the identity ma- 
trix. But A = [ ^J* ~ ) , which is clearly not skew- 
Hermitian. 

If two matrices are unitarily equivalent, then the first 
is Hermitian if and only if the second is Hermitian. Thus 
a sufhcient condition for the matrix A = V XTd\a,g{—(l)j + 
2iTkj)Xj^^V^ to be Hermitian is for the matrix Xt to be 
unitary, so that the Hermitian matrix diag(— 0j -1- 2Tikj) 
would be unitarily equivalent to A via the unitary V Xt- 
This condition is not in general necessary. The neces- 
sary and sufficient condition is that X^rpXT commutes 
with diag(— -|-27rfcj), which can easily be derived, once 
more using the stated matrix property. Hence, the inter- 
section in Equation [lO] can be removed if the statement 
[X^rj^Xx^ diag(— 0j -f 2T:kj)] = is appended to the condi- 
tions. 

We may count that there are four freedoms at our dis- 
posal when picking a member of Q-u . The first comes from 
our characterization of Qu and is the choice of any block- 
diagonal unitary U in the group U{pf , . . . ,Pmo)^ whose 
dimension as a manifold is determined by the amount of 
degeneracy in the eigenvalues of the state p- (or equiva- 
lently, This dimension can range from (d— 1)^, when 
P- is a pure state, to d, when no two of p-'s eigenvalues 
are equal. 

Once U has been chosen, the second freedom lies in the 
choice of diagonalization, {V,T), of J7, which is charac- 
terized by the set Vu- As we have seen, choosing {V,T) 
comes down to choosing a unitary in U{p^ , . . . ), as 
well as any permutation matrix P. Thus the size of Vjj 
is determined by the amount of degeneracy in the eigen- 
values of U. 

The third freedom comes in the form of d choices of 
integers for ki through k^. 

The final freedom is in the choice of member from the 
centralizer group Cgl{T). The size of this freedom is 
dependant on the number of unique values contained in 
T, i.e. the degeneracy of U once again. We must also 
remember that this freedom is restricted if we wish for 
our Hamiltonians to be Hermitian. 

We now turn toward an illuminating example. Recall 
that Cc is defined to be the Lie algebra generated by 
our system controls. In the case where Q-u and iLc have 
empty intersection, we simply do not have enough control 
to do as we wish On the other hand, if Q-u C iCc, 
then we have too much control for our scheme to be use- 
ful; we can evolve under Hd directly. And so for us, the 
only interesting case is when and iCc intersect, but 
neither set is contained within the other. In this case, 
perhaps, Hd lies outside of iCc, but we can effectively 



simulate the piece we are interested in for a specified 
time ^0 with some H G iCc H Q-h . 

It is not immediately clear that this third and most 
interesting case should ever happen. Perhaps each H 
in Q-j-i is implementable if and only if one H in Q-u is 
implcmcntablc. The following provides a counterexample 
to this idea. 

Suppose we are interested in passing information from 
one end of a four qubit spin chain to the other, and 
that this spin chain's internal Hamiltonian is given by 
a dipole-dipole interaction 

^ J,,(4'=)a(')+af)4')-2aWa«) (11) 

l<fc<i<4 

where Jki oc \k — l]^"^ falls off as distance cubed. 

Suppose also that we have some control over the sys- 
tem, the global X and Y rotations whose Hamiltonians 
are given by 

4 4 

= E ^-''^ and Hy = 4'^ (12) 

k=l k=l 

The propagation of quantum information along spin 
chains in the case of either no or limited control has 
been well studied [H]. The simplest method, which nat- 
urally yields poor fidelities, is to allow the free evolution 
of an XY nearest-neighbour Hamiltonian with uniform 
couplings and wait for a specific amount of time |14) . 
If one has the ability to engineer non-uniform nearest- 
neighbour couplings, then the same idea can produce unit 
fidelity transfers [13 US] . Again in the case of no control 
on the spin chain, if one has a dual rail spin chain, or 
equivalently, a single rail of qutrits, unit transfer fidelity 
of a qubit can be achieved with a wide class of chain 
Hamiltonians [T7| . If one adds the ability to perform cer- 
tain controls to a nearest-neighbour chain, then pulse se- 
quences can be found on uniformly coupled models with 
unit fidelity transfers [18] , and with engineered couplings 
the state of the chain doesn't even need to be initial- 
ized [12] • The example presented here excels in that the 
chain's interaction is not restricted to nearest-neighbour 
couplings, and that the controls are only required to be 
global rotations, but lags in its numerical tractability for 
long chains. 

We take the following XY Hamiltonian as our desired 
Hamiltonian 

HxY = X: C,iai'^alt+'^ + Vf+i)) (13) 

k=l 

because if the coefficients are chosen as Ck — 
^■y/A:(4 — fc), then information can be passed from the 
first spin on the chain to the last spin on the chain with 
unit fidelity simply by waiting for a time 7r/A|15|. In 
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other words, for any a and /? we have the following prop- 
erty: 



a\nu)+P\um 
= a luit) + P \mi) 



(14) 



By taking nested commutators of our three control 
Hamiltonians {Hj., Hy, and Hjjjj) until all further nested 
commutators lie in the span of the previous ones, we find 
that our Lie algebra has a dimension of 96, with the di- 
mension of the full control space being 4^ — 256 [9 . The 
norm distance between Hxy and its projection onto the 
space iCc can be computed to be non-zero, so that Hxy 
is not in our control algebra. This means that if a so- 
lution is found to our problem, it will be a non-trivial 
one. 

Let our initial state be the pure state pi = 
Itm) (tmL and our fraction of internal evolution be 
77 — 0.95. We choose the diagonalization unitary [/_ so 
that D = diag(l,0, • • • ,0), hence Gu = U+U{1, 15)C/1. 

We wish for our solution to retain the property in 
Equation |14[ but in choosing a unitary U from Qu , we 
can only be guaranteed that this property holds for one 
state, namely, UintUUint Itiii) = Uiit)- To overcome 
this limitation, we use the following trick. We choose 
a subspace of Cc, denoted Ccii), such that each ma- 
trix in this subspace has as an eigenvector. We 
can expect such a subspace to exist since IJ-iiJ,) is an 
eigenvector of both Hint and [H^jHy]. In our calcula- 
tions, such a subspace was found having a dimension of 
9, where no attempt was made to find the largest such 
subspace. Now if we were to successfully perform the 
optimization 



H. 



calc 



argmin 1 1 _ff 



{H)\\ 



(15) 



where iVcdl) is the projection operator onto iCc{i), we 
would be ensured that 



Uinte 



^j(l-J))t0ffca!c 



u,nt \mi) = \mi) (16) 



and hence by linearity, and since Hcaic G Gu^ the desired 
property 

f/mte-^(i-'')*''^--C/,„t(a im;) + 13 mi)) 

= a |;Ut> + 13 mi) (17) 

would be recovered for any a and (3. 

Actually performing this minimization is non-trivial 
since the spaces we must optimize over are non-linear, 
and moreover, a full optimization would require a double 
optimization, first over Qu and then over Tyjj, and 
CcLiT). However, in our example a partial optimization 
turns out to suffice, where we optimize only over Gu, 
leaving our other freedoms fixed at simple values (kj — 
for all j, Xt = and letting Matlab's diagonalization 
function choose the member of 'Djj)- 




I exp ( -itHxY ) I TJ-J-i) 

■ (J-J-iT|exp(-itHi„t) iTii-l) 

■ (J"U.T|exp(-ilJj Ho(T)dT) |TJ,ii> 

|exp(-il|, Ho(x)d-C) |iJ,ii) 



FIG. 2: The fidelity of information transfer from the first 
to the last qubit on the spin chain for the states of interest, 
Iti^^) a-nd |4,i4-4')i evolving under the numerically discovered 
Hamiltonian Hq. Also included for comparison are the fideli- 
ties of transfer under the Hxy and Hint Hamiltonians. Recall 
that 14^44) is an eigenvector of all Hamiltonians involved. 



This scheme yields a Hamiltonian Hcaic which lies in 
both Qi-i and iCc{\-), and obtains the global minimum 
of for the cost function written down Equation 15 It 
should be noted that both of the maximizations 



H'^^i^ = argmax | (lUtI e-'^i-'')*"^ mU) | 
HeCcii) 

H'J,i, = argmaxd (U;tl e-'^'-^'^*"" lUU) I 
HeCc 

+|(;;;;|e-»(i-")*««i;u;)|) 



(18) 
(19) 
(20) 



would be much simpler, and produce effectively the same 
result in this example. The difference, however, is that 
these two methods demand that our operation be imple- 
mentable while seeking a maximum in state transfer fi- 



delity, whereas the minimization in Equation 15 demands 
that the fidelity be perfect while attempting to minimize 
the distance to implementability. 

Figure [2] shows the fidelity of state transfer as a func- 
tion of time, where Hq is the piecewise constant Hamil- 
tonian which begins and ends as Hint, and is equal to 
Hcaic for a duration (1 — 77)^0 in between. Similar results 
have been computed for both 5 and 6 qubits. 

Having calculated the Hamiltonian Hcaic, the methods 
of optimal control theory can be used to generate control 
sequences which implement it [2D]. If desired, the con- 
trol sequence can subsequently be concatenated with its 
time and phase reversal to generate a universal rotation 
gate pi] . 

In conclusion, we have provided a complete characteri- 
zation of the set of Hamiltionians which simulate the evo- 
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lution of a given density matrix under a desired, but not 
implementable, Hamiltonian. Moreover, the above ex- 
ample demonstrates that in a non-trivial physically moti- 
vated scenario, it is possible for this set to have nonempty 
intersection with the control algebra, thereby providing a 
realizable means of doing what can not be done directly. 
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